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ABSTRACT 

We  are  exploring  the  use  of  acoustic  resonant  techniques  to  identify  defective 
components  by  monitoring  the  fundamental  modes  of  vibration  that  arise  following  an  impulse 
load.  We  employ  the  spectrum  of  the  filtered  data  as  a  feature  for  use  in  automated  statistical 
pattern  recognition.  The  feature  provides  a  measure  of  the  severity  of  the  defect  using  the 
acoustic  signature  of  a  defect  free  component  as  a  baseline.  Defects  change  the  signature  by 
shifting  or  splitting  the  frequencies  and  introducing  beats.  Identifying  these  changes  in  the 
spectrum  relative  to  baseline  data  is  difficult  because  of  the  large  number  of  natural  frequencies 
associated  with  complex  geometries  and  the  presence  of  noise  in  the  data.  In  the  case  of  real¬ 
time  health  monitoring,  obtaining  a  representative  spectrum  is  challenging  since  averaging 
cannot  be  employed  without  the  use  of  multiple  transducers.  We  have  developed  a  new  approach 
to  obtain  a  smoother  spectrum  by  analyzing  sequential  data  using  an  approach  similar  to 
Rayleigh’s  quotient  method  for  eigenvalue  problems  and  super-resolution  techniques.  The 
resulting  spectrum  enhances  the  natural  modes  of  vibration  and  provides  a  better  feature  for 
automated  classification  techniques.  The  numerical  implementation  of  the  algorithm  requires 
modest  computational  resources.  The  algorithm  for  generating  the  Rayleigh  quotient  has  been 
successfully  tested  on  numerically  generated  waveforms  and  experimentally  acquired  data.  In  all 
cases,  there  was  a  significant  reduction  in  the  inherent  noise  and  the  numerically  generated  noise 
from  digital  filtering. 
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BACKGROUND 

The  derivation  of  Rayleigh’s  quotient  is  outlined  below  to  provide  insight  into  how  it  can  be 
obtained  as  sequential  transforms.  These  transforms  can  then  be  used  to  simplify  the  numerical 
techniques  so  the  approach  can  be  implemented  using  modest  computational  resources. 

Consider  a  signal  G- , 


G,  =  A  exp(-2^/0/t;  )  =  A  exp(-2^/0A//),  i  =  0,  N  - 1  (1) 

where  A{  is  the  amplitude,  /0  is  a  single  frequency,  and  Ms  a  time  vector.  The  Discrete  Fourier 
Transform  (DFT)  of  G.  is: 


i  i=N- 1  7  1  N- 1 

Ht=  —  y  G.  exp(2^i— )  =  —  y 
1  N  U  N}  NU 

2  7ll 

Let  a  =  -2 nf0At,  P  =  — — ,  a  =  exp (ja). 
Let  Ai=  1,  and  sum  the  series  to  give  //,  as: 


A  exp(2^7(T  -  f0At)),  l  =  0,..,N- 1 
b=txV(jp)  b=b(i),  e=^- 


(2) 


\n  sin(#)  j  (3) 

Equations  (2,3)  are  the  usual  Fourier  transform  of  G- .  Equation  (3)  can  also  be  obtained  by  constructing 

a  covariance  matrix  using  super-resolution  methods  [4].  A  set  of  vectors  is  defined,  using  a  step-wise 
approach,  as: 


/  N- 1  N- 2  n 

xl={a  a  •••!) 
x2  =  ( aN~ 2  aN~ 3  •  •  •  0),  etc 

x^  =  (1 0  0 . 0) 

B  =  (bN  l  bN~2---l)T 


The  Rayleigh’s  quotient,  XN ,  is  [5]: 

XN  =  (XN  XN- 1  XV-2 . Xl) 

xnB  =  (xnB  xn_xB  xn_2B  '-xxB)t  (5) 

_  ( xnB)t(xnB ) 

* 

Cov(xN )  =  (x^)  (x^y  )  (6) 
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Equation  (6)  is  the  covariance  matrix  of  source  vector  and  B  is  the  test  vector.  Rayleigh’s  quotient 
provides  an  alternate  to  the  Fast  Fourier  Transform  (FFT),  but  requires  a  large  matrix  and  is  numerically 
intensive.  Figure  1  shows  MATFAB  [5]  code  used  to  generate  the  quotient.  The  complexity  can  be 
reduced  by  analyzing  the  structure  of  the  quotient. 

function  RAFY  l=Rayleigh_old(sig,t,freq_F) 
vector  l=sig; 

NN=length(sig); 
xx  1  =zeros(NN,NN); 
for  lag=0:NN-l 

xxl(:,lag+l)=[zeros(l,lag),  conj(vectorl(l:NN-lag))]f; 
end 

R1  =  xxl’*xxl;  %  Correlation  matrix 

fO=freq_F; 

for  k=l:NN 

S=exp(-2*pi*i*f0(k)*t); 

RAFY  1  (k)=(S*Rl  *S’)/(S*Sf); 

end 

return 

Figure  1  MATLAB  code  for  generating  the  Rayleigh  Quotient 


It  should  be  noted  that  in  above  definition  XN  has  the  structure  of  an  eigenvalue  but  it  is  actually  a 
transform  of  the  original  signal.  It  has  an  implicit  transform  variable  /  =  0,..,  N  —  1 .  The  following  shows 
how  XN  can  be  given  as  the  sum  of  the  squares  of  sequential  Fourier  transform. 


Consider  only  one  term  in  (5): 


K  = 


BtB 


1 


a\N- 1  .  sa\N- 


a\N- 1  .  sa\N- 


BLB  \  b 


(-f-1 + (-)*- +•••+!  +  {-y-1  h — hi 


1  1 

J 

1 

V  J 

\aN-bN \ 

N  aN~l  | 

a-b 

k  J 

bN -1 1 

1 

i  i 

N  ( ab)N~l 


a^-b^ 
a-b 


1  Sill  2N9  ,,|rrl2 

- 5 —  =  N\H,\ 

N  sin2  9 


(7) 


Similarly, 

_  (x2 Bfjx^B ) 

D^D 


(8) 
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Considering  the  sum  of  all  the  terms  and  the  identity  (< a  —  b)T  (a  —  b)  =  4  sin2  0 : 


•N- 1 


K  = 


N(a-b)T 

1  1 
iV  4 sin2# 


1  1 


f  ,  >\ 

a-  b 

a2  -b2 


yaN -bN  j 


bN (a-b) 

{4sin2#  +  4sin22#  +  4sin23#  + — h  4  sin2  TV#} 


(9) 


Using  the  summation  formula  we  have: 

(10) 

sin#  J  47V sin2  #[  (27V  +  l)sin#j 

Equation  10  is  the  Rayleigh  Quotient  transform  of  a  signal  with  a  single  frequency  fO  as  compared  to 
equation  2  and  equation  3  which  are  Fourier  transforms.  For  comparison  equation  2  has  been  simplified 
for^4,=l  as: 


SinN9_ 

TV  sin# 


(11) 


Equations  (10)  and  (11)  incorporate  sine  functions  and  for  the  ranges  involved,  can  be  less  than  1. 
Therefore,  the  Fourier  transforms  can  be  zero  resulting  in  negative  spikes.  The  Rayleigh  quotients  will 
never  be  zero. 


For  numerical  implementation,  the  Rayleigh  quotient  can  be  represented  as  the  sum: 

/L  =  — - h —  (4  sin2  #  +  4  sin2  2#  +  4  sin2  3#  H —  +  4  sin2  TVi 

N  TV  4sin2  # 

=  —{h2in  +  H2,n-i---H2,o)} 

TV  1  J 

=  J.//r(.v:)  ,  ,//r(.v2)  #2(xw)l 

[  N2  N2  N2 


where  H f(N  —  /)| 


Sin(N^  -  /)# 

TVsin# 


i  =  1,2,3-  •• 


(12) 
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Figure  2  is  the  corresponding  MATLAB  code  for  generating  the  sequential  transform  which  uses  the  FFT 
routine.  The  input  is  the  time  series  signal  and  the  output  is  Rayleigh  quotient. 


function  Ys=Rayleigh(x) 

tot=0; 

x=fliplr(x); 

for  i=0:length(x)-l,tot=tot+abs(fft([x(i+l:end),zeros(l,i)])).A2;end 

Ys=tot/length(x); 

return 

Figure  2  code  for  efficient  generation  of  the  Rayleigh  Quotient. 
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RESULTS 

The  algorithm  for  generating  the  Rayleigh  quotient  was  tested  on  numerically  generated 
waveforms  and  experimentally  acquired  data.  In  all  cases,  the  inherent  noise  and  the  noise  generated 
numerically  (i.e.  sidelobes  from  numerical  filters)  was  reduced  significantly.  Figure  3  shows  the 
normalized  FFT  of  a  signal  and  figure  4  shows  the  transform  obtained  with  Raleigh  quotient.  The  signal, 
0.5sin(2;zT000T)  +  1.0sin(2;rl500T) ,  was  generated  numerically  and  corrupted  with  random  numbers 
with  a  signal  to  noise  ratio  of  -14  dB. 


1 

0.8 
-g  0.6 

'c 

CD 


E  0.4 


frequency  (Hz) 

Figure  3.  fft  of  0.5sin(2;rl0000  + 1 .0sin(2;rl 5000 


Figure  4.  Rayleigh  Quotient  of  0.5sin(27rl000f)  +  1 .0sin(27rl500^) 


We  are  currently  employing  a  compact  tensile  specimen  in  a  systematic  study  to  quantify  the 
propagation  of  a  crack  in  terms  of  changes  in  the  fundamental  modes  of  vibration.  A  crack  was  induced 
in  a  tensile  specimen  and  grown  in  regular  increments.  Raw  vibration  data  was  collected  at  each 
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increment  using  two  accelerometers  mounted  to  the  specimen,  an  instrumented  hammer  to  generate  the 
impulse  load,  and  a  spectrum  analyzer  to  collect  the  data.  Approximately  100  data  sets  from  each 
accelerometer  were  collected  for  each  specimen.  Figure  5  shows  the  average  of  the  FFT  for  second 
fundamental  (9407  Hz,  defect  free)  at  each  test  increment.  Results  are  similar  for  the  first  and  third 
modes.  Figure  6  shows  the  second  fundamental  computed  for  a  single  data  set  using  the  FFT  and  method 
of  sequential  transforms.  The  data  was  filtered  using  a  lowpass  filter  with  a  cutoff  frequency  of  12  kHz. 
The  figure  shows  how  the  method  of  sequential  transforms  removes  sidelobes  and  sharpens  the  center 
frequency. 


Figure  5.  Second  fundamental  frequency 
for  a  compact  tensile  specimen  with  different  crack  lengths. 


Figure  6.  FFT  and  Rayleigh  Quotient  for  measurements  from 
two  accelerometers  mounted  on  the  compact  tensile  specimen. 
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SUMMARY  AND  CONCLUSION 


We  have  introduced  a  new  transform  method  for  obtaining  a  smoother  spectrum  by  analyzing  sequential 
data  using  an  approach  similar  to  Rayleigh’s  Quotient  method  for  eigenvalue  problems  and  super¬ 
resolution  techniques.  This  provides  a  means  for  improving  the  quality  of  features  for  automated 
classification  techniques  without  the  need  for  acquiring  multiple  samples  or  using  additional  transducers. 
We  have  derived  the  approach  analytically  and  provided  a  numerically  efficient  means  to  obtain  the 
transform.  The  algorithm  has  been  successfully  tested  on  numerical  data  with  artificially  induced  noise 
and  on  experimental  accelerometer  measurements. 
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